LPT-Orsay-10-22 



One-loop contribution to the neutrino mass matrix in NMSSM with right-handed 

neutrinos and tri-bimaximal mixing 

Debottam Das* 
Laboratoire de Physique Theorique, UMR 8627, 
Universite de Paris-Sud 11, Bailment 210, 
91405 Orsay Cedex, France 

Sourov Royt 

Department of Theoretical Physics and Centre for Theoretical Sciences, 
Indian Association for the Cultivation of Science, 
2A & 2B Raja S.C Mullick Road, Kolkata 700 032, India 
(Dated: August 5, 2010) 

Neutrino mass patterns and mixing have been studied in the context of next-to-minimal super- 
symmetric standard model (NMSSM) with three gauge singlet neutrino superfields. We consider 
the case with the assumption of R-parity conservation. The vacuum expectation value of the singlet 
scalar field S of NMSSM induces the Majorana masses for the right-handed neutrinos as well as 
the usual /i-term. The contributions to the light neutrino mass matrix at the tree level as well as 
one-loop level are considered, consistent with the tri-bimaximal pattern of neutrino mixing. Light 
neutrino masses arise at the tree level through a TeV scale seesaw mechanism involving the right- 
handed neutrinos. Although all the three light neutrinos acquire non-zero masses at the tree-level, 
we show that the one-loop contributions can be comparable in size under certain conditions. Pos- 
sible signatures to probe this model at the LHC and its distinguishing features compared to other 
models of neutrino mass generation are briefly discussed. 



PACS numbers: 12.60.Jv, 14.60.Pq, 14.60.St, 
I. INTRODUCTION 

Several mechanisms of the generation of neutrino 
masses and mixing in the context of a supcrsymmetric 
model have been explored in various works. One of the 
most popular attempts in this direction is to relax the as- 
sumption of R-parity conservation in the minimal super- 
symmetric standard model (MSSM) by including explicit 
bilinear and/or trilinear R-parity violating interactions in 
the superpotential and the scalar potential [1, 2]. One can 
also consider models with spontaneous R-parity violation 
[3-5] via a singlet sneutrino vacuum expectation value. 
The low energy limit of such models, where the singlet 
sneutrino field is decoupled, can be thought of as the bi- 
linear R-parity violating scenario. Thus there are several 
possibilities within the context of R-parity violation in 
MSSM. In fact, each of them has been studied in detail 
in connection with the observed neutrino mass patterns 
and mixing as provided by the neutrino oscillation exper- 
iments. The possible collider signatures of R-parity vio- 
lating models have also been studied in great details and 
correlation between neutrino mixing angles and the decay 
branching ratios of the lightest supersymmetric particle 
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(LSP) have been obtained [6-15]. 

Another interesting and well studied procedure of 
small neutrino mass generation in a supcrsymmetric 
model, with the observed mixing pattern, is the seesaw 
mechanism[16, 17] with the introduction of right-handed 
neutrino superfields [18-20]. In order to generate small 
neutrino masses, one introduces AL = 2 heavy Majo- 
rana mass terms in the superpotential in addition to the 
trilinear lepton-number conserving Yukawa interactions 
involving the right-handed neutrino superfields. As long 
as the neutrino Yukawa couplings are of order one, light 
neutrino masses ~ 10 -2 eV require the Majorana masses 
to be ~ 10 15 GeV or so. However, such a high seesaw 
scale is difficult to probe at the LHC or future linear 
collider experiments. A viable alternative is to look at 
Te V -scale seesaw mechanism where small active neutrino 
masses are generated with the help of neutrino Yukawa 
couplings as small as 10 -6 (same as the electron Yukawa 
coupling) and this makes the Majorana mass scale of the 
right-handed neutrino of the order of ~ TeV plausible. 
This gives one an opportunity to test the seesaw models 
at the LHC. The signatures of TeV scale supersymmetric 
seesaw models will be briefly outlined later along with a 
discussion of the signatures of R-parity-violating models. 

On the other hand, MSSM is plagued by the so-called 
"/^-problem" which asks the question that why the scale 
of the supersymmctry preserving /z-term should be of the 
same order as the soft supersymmetry breaking terms, 
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which are of the order of TeV. One of the possible so- 
lutions to this problem is the next-to-minimal super- 
symmetric standard model (NMSSM) , where a standard 
model singlet supcrfield (S) is introduced to the MSSM 
superfields with a coupling XSH u Hd in the supcrpoten- 
tial (for review and phenomenology see [21, 22]). The 
scalar component of S gets, in general, a non-zero vac- 
uum expectation value (VEV) of the order of ~ TeV, 
as long as the soft mass parameters corresponding to 
the singlet scalar field arc in the same range. This 
solves the "//-problem" because the //-parameter gen- 
erated in this way has the right order of magnitude if 
one considers a coupling A ~ 0(1). In order to gen- 
erate active neutrino masses and appropriate mixing in 
the neutrino sector one either includes R-parity viola- 
tion in the superpotcntial[23, 24] and the scalar poten- 
tial or introduces gauge-singlet neutrino superfields Ni 
with appropriate couplings with the MSSM superfields 
and the singlet supcrfield £[25]. In the latter case, the 
gauge-singlet neutrino superfields Ni can have Majorana 
masses around the TeV scale if there is a coupling of the 
type nNfS in the superpotcntial. When the scalar com- 
ponent of S gets a VEV of the order of TeV scale, the 
right handed neutrinos also acquire an effective Majorana 
mass around the TeV values as long as the dimcnsionlcss 
coupling k is order one [25]. Here it is assumed that the 
superpotential has a discrete Z3 symmetry which forbids 
the appearance of bilinear terms in the superpotential 
[26]. 

In this study, within the framework of this TeV scale 
seesaw model mentioned above, we calculate the one-loop 
contributions to the neutrino mass matrix with R-parity 
conservation and study the effect of these contributions 
to the neutrino mass patterns and mixing angles. In 
other words, we consider the case where only the scalar 
field corresponding to the singlet superfield S gets a non- 
zero VEV along with the neutral Higgs fields. We will 
show later that these one-loop contributions can be sig- 
nificant and can change the region of parameter space 
allowed by the three-flavor global neutrino data in com- 
parison to the tree level results. 

The plan of the paper is as follows. In Sec. II we will 
provide a discussion on the three-flavor neutrino mixing 
and illustrate the general pattern of our analysis that we 
are going to follow. Sec. Ill describes the model along 
with the minimization conditions of the neutral scalar 
potential. One-loop contributions to the neutrino mass 
matrix in the R-parity conserving scenario and the re- 
sulting neutrino mass patterns, which satisfy the three 
flavor global neutrino data, are discussed in Sec. IV with 
numerical results. In Sec.V we outline the possible ways 
to probe this model at the LHC and present a short crit- 
ical discussion of the signatures of neutrino mass models 
involving spontaneous and/or bilinear R-parity violation. 
We summarize in Sec. VI with possible future directions. 



II. NEUTRINO MIXING 

The solar, atmospheric, accelerator, and reactor neu- 
trino experiments have shown strong evidence in favor of 
non-zero neutrino masses and mixing angles [2 7]. In addi- 
tion, there is an upper bound on the sum of neutrino mass 
eigenvalues ~1 eV from cosmological observations [28]. 
The bound on the 11-element of the neutrino mass ma- 
trix resulting from the non-observation of neutrinoless 
double beta decay is <0.3 eV[29]. The global 3-flavor fits 
of various neutrino oscillation experiments point toward 
the following 3ct ranges of the neutrino oscillation param- 
eters, namely the two mass-squared differences and three 
mixing angles [30]: 

sin 2 0i2 = 0.25 - 0.37, sin 2 23 = 0.36 - 0.67, 
sin 2 6>i3 < 0.056 

Am|i = (7.05 - 8.34) x 10" 5 eV 2 , 

|Am^| = (2.07- 2.75) x 10" 3 eV 2 , (1) 

where Am?- = to 2 — to 2 . One can see from these numbers 
that there are two large mixing angles and one small 
mixing angle among the three light neutrinos with a mild 
hierarchy between the mass eigenvalues. 

The three flavor neutrino mixing matrix U can be 
parametrized as follows, provided that the charged lep- 
ton mass matrix is already in the diagonal form and the 
Dirac as well as Majorana phases are neglected: 



(C12C13 S12C13 S13 ) 

-S12C23 - C12S23S13 C12C23 - S12S23S13 S23C13 
S12S23 - C12C23S13 -C12S23 - S12C23S13 C23C13 J 

(2) 

where cy = cos0ij, = sin0,j and i,j run from 1 to 3. 

The mixing angle data coming from solar, atmospheric 
and reactor sector indicate that 0i2 ~ 34°, 02 3 « 45°, and 
$13 < 13°. This is popularly known as the bilarge pattern 
of neutrino mixing. In order to understand the conse- 
quences of such mixing in the zeroth order, one considers 
the tri-bimaximal structure of the neutrino mixing[31] 
where 2 3 = f , $13 = and sin0i 2 = 4j. 

With this tri-bimaximal pattern, the unitary neutrino 
mixing matrix turns out to be 
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Considering mi, 11J2 and m 3 as the three light neutrino 
mass eigenvalues, we use the matrix U v to obtain the 
neutrino Majorana mass matrix in the flavor basis as 
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\{—mi + m 2 ) 
g(mi + 2m 2 + 3m 3 ) 
■(—mi — 2m 2 + 3m 3 ) 



|(mi - to 2 ) 
g(— mi — 2m 2 + 3m 3 ) 
+ 2m 2 + 3m 3 ) 



(4) 



r 



We can see that a particular structure of neutrino mass 
matrix emerges from the requirement of tri-bimaximal 
mixing, in terms of the neutrino mass eigenvalues. Given 
a specific model for generating the neutrino mass ma- 
trix, one can easily connect the model parameters with 
the neutrino mass eigenvalues with the help of Eq.(4). 
This way one can study the normal, inverted or quasi- 
degenerate mass pattern of the light neutrino mass eigen- 
values and try to see the requirement on the model pa- 
rameters to produce the tri-bimaximal pattern of neu- 
trino mixing. In this work, we will try to explore 
the next-to-minimal supersymmetric standard model 
(NMSSM) where neutrino mass is generated because of 
the introduction of three right-handed neutrino super- 
fields with the possible interaction terms. Though the 
assumption of tri-bimaximal mixing in the neutrino sec- 
tor is not generic, in the present context it is quite illus- 
trative in studying the role of the soft SUSY breaking pa- 
rameters on the neutrino mass eigenvalues. At the same 
time, the acceptable domain of the soft parameters con- 
sistent with neutrino mass eigenvalues and tri-bimaximal 
mixing angles would hardly change with any small shift 
in 0i 3 . 

As mentioned in the introduction, this model was pro- 
posed in Ref. [25] where the case with spontaneous vio- 
lation of R-parity was studied with possible implications 
on neutrino mass eigenvalues and mixing angles at the 
tree level. In the present study we shall consider the case 
when R-parity is conserved and the neutrino mass gener- 
ation at the tree level is entirely due to the seesaw mech- 
anism involving the TeV scale right handed neutrinos. 
Our aim would be to see if this model can produce the 
acceptable neutrino mass eigenvalues and mixing angles 
when the neutrino mass matrix receives contributions at 
the tree as well as one-loop level. An attractive feature of 
this model is that, the right handed sncutrino in the form 
of LSP may become a valid cold dark matter candidate 
of the universe [32]. 

This model can also accommodate spontaneous CP 
and R-parity violation simultaneously. In that case, the 
neutrino sector is CP violating and the resulting effects 
on the neutrino masses and mixing angles were studied in 
Ref. [33] . Similarly, spontaneous R-parity violation moti- 
vated by a flavor symmetry may produce tri-bimaximal 
mixing pattern in the neutrino sector [34]. However, in 
the present context we consider the case where neutrino 
sector conserves CP symmetry along with R-parity. 



There have been some other studies which address the 
neutrino experimental data in some other extensions of 
NMSSM. One of these proposals is discussed in Ref. [23], 
where the effective bilinear R-parity breaking terms are 
generated through the vacuum expectation value of the 
scalar component of the singlet superfield S. In this 
case, only one neutrino mass is generated at the tree 
level whereas the other two masses are generated at the 
one-loop level. In another model[24], non-zero masses for 
two neutrinos are generated at the tree level by including 
explicit bilinear R-parity violating terms along with the 
R-parity breaking term involving S. It is interesting to 
note that, the R-parity violating NMSSM model may of- 
fer a valid dark matter candidate in the form of gravitino 
as the R-parity violating decay channels of the gravitino 
are extremely suppressed because of weak gravitational 
strength[35]. 



In another class of models, gauge-singlet neutrino su- 
perfields were introduced to solve the //-problem, which 
can simultaneously address the desired pattern of neu- 
trino masses and mixing [36]. The detailed study of neu- 
trino masses and mixing in this model was presented in 
Ref. [37] and the correlations of the lightest neutralino de- 
cays with neutrino mixing angles were discussed. Subse- 
quently the dominant one-loop contributions towards the 
tree level neutrino masses have also been presented [38]. 
Similar analyses for one and two generations of gauge- 
singlet neutrinos were presented in Ref. [39] and some 
other phcnomcnological implications, in particular the 
possible signatures at LHC were addressed. Neutrino 
masses consistent with different hierarchical scenarios 
and tri-bimaximal neutrino mixing can also be gener- 
ated in an R-parity violating supersymmetric theory with 
TeV scale gauge singlet neutrino superfields, where the 
^-term was not generated by the vacuum expectation 
values of the singlet sneutrino fields [40]. Another in- 
teresting avenue in this direction is to study the role of 
possible higher dimensional supersymmetry breaking op- 
erators in the hidden sector which may render the TeV 
scale soft SUSY breaking trilinear and bilinear couplings 
involving the sneutrinos to produce the observable mass 
and mixing angles for the neutrinos[41]. 
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III. THE MODEL AND MINIMIZATION 
CONDITIONS 

In this section we review the model along the lines of 
Ref.[25] and discuss its important characteristics. We 
introduce the singlet superfield S along with three right- 
handed neutrino superficlds Ni. The superfields Ni 
are odd and the superfield S is even under R-parity. 
The most general supcrpotcntial consistent with R-parity 
conservation is 



(5) 



where 
Wnmssm 



W = W-NMSSM + Wsmglet , 

= ft(H d Q l )b l +mQ l H u )U J 



f?(H d Li)Ei + \ H (H d H u )S + -^S 3 , (6) 



Wsinglet = f^LiA^Nj + ^N?S . 



(7) 



Here H d and H u are down-type and up-type Higgs super- 
fields, respectively. The Qi are doublet quark superficlds, 
Uj[Dj] arc singlet up-type [down-type] quark superfields. 
The Li are the doublet lepton superficlds, and the Ej 
are the singlet charged lepton superfields. The indices 
i,j = 1,2,3 are generation indices. Note that we have 
imposed a Z3 symmetry under which all the superficlds 
have the same charge. This symmetry forbids the ap- 
pearance of the usual bilinear ^-tcrm in the supcrpotcn- 
tial. The /Lt-term is generated spontaneously through the 
vacuum expectation value of the singlet scalar S. In a 
similar way soft supersymmetry breaking potential can 
be written as 



soft 



t/NMSSM 
^soft 



^Single 



where V r s X SSM includes the MSSM soft supersymme- 
try breaking terms along with a few additional terms as 
shown below: 



t/NMSSM 
''soft 



V s ^ M + mi\Sf + 
A, 



A X H H d H u S + Am-^S 6 + H.c 



(9) 



The term Vsingiet is composed of the soft masses and the 
trilinear interactions corresponding to the fields Ni - . 



Vsingiet = ™%n. l^| 2 + 



A v fy j L i H u N j + Am'-^SNf + H.c 



(10) 



We have taken a common trilinear coupling A for the 
singlet fields N and S and to is a mass scale. In a su- 
pcrgravity motivated scenario, it is a common practice 
to choose to = TO5 = mfjf{* and also a universal trilin- 
ear parameter for the fields 5", N. Since these fields are 
gauge singlet, we assume such universality to hold also 
at the electroweak scale. Similarly, the mass parameters 
to 5 and nij-jjy, are very much insensitive to Rcnormal- 
ization Group Equation (RGE) running and their values 
at the weak scale can be taken to be the same as the 
values at the high scale. In addition, we have chosen all 
the parameters ff, ff, Ajvj, A#, A s ,/™- and f"- to be real. 

The scalar potential of this model can be written as 



V = V F + V D + V soit 



(11) 



(8) where the neutral part of Vf and Vd can be written as 



V, 



neutral 



V "neutral 
n 



i 

+ J2\f^iH°) + X Nj N j S 



i 

r 



KDiNj + \ H H° d S 



\ H (Hm + ^N? + ^-S 2 



(12) 
(13) 



In the above, the repeated indices always mean to sum 
over the generations. However, the summation sign is 
used in special cases if required. The VEVs are deter- 
mined by the minimization of the potential (vide Eq.(ll), 
(12) and (13)). Here we explore the possibility when 
only scalar component of the gauge singlet superfield 



S acquires a VEV along with the doublet Higgs fields. 
The right-chiral sncutrino N can only have a vanishing 
VEV and thus R-parity is unbroken. On the other hand, 
when the right-chiral sneutrino N acquires a VEV then 
R-parity is spontaneously broken and an effective bilinear 
R-parity violating term of the form 6iLiH u is generated, 
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where = f v (Ni). However, the case of spontaneous 
R-parity violation will be studied in a separate work [42]. 
Note that, a global continuous symmetry such as lepton 
number cannot be assigned to the supcrpotcntial involv- 
ing the singlets S and N% . Thus this model is completely 



dV 
dv d 

+ 

dV 

+ 
+ 

dV 

dv s 

+ 

dV 
dV 

9v Ni 

+ 



Here g\ and gi are the U(l) and SU(2) gauge couplings, 
respectively, and tan/3 = v u /v d . fhi is the soft SUSY 
breaking mass parameter of the left chiral sneutrinos. We 
have assumed that the neutral scalar fields can develop, 
in general, the following vacuum expectation values 

v d = (H° d ); v u = (H° u ) ] v s = (S); 
vs, = {i/i); v$. = (Ni). (15) 

As has already been mentioned, in the present context 
we will consider the solutions va. = vp t =0 and v s =/= 
to analyze the neutrino spectra. In our subsequent dis- 
cussion, we will also ignore the terms in the minimization 
equations which are bilinear in the neutrino Yukawa cou- 
plings. Note that in order to generate very small masses 
for the active neutrinos ( $ 0.1 eV) using this TeV scale 
seesaw mechanism, the neutrino Yukawa couplings (/") 
should be below 0(1O -6 ), which is around the magnitude 
of the electron Yukawa coupling. 

The VEV v s comes out as the solution of the following 
cubic equation (neglecting the Yukawa term), 

\ 2 s Vs + Am\ s v 2 s + 2v s (m 2 s + X 2 H v 2 + X 2 H v d + X H \ s v d v u 
+ X 2 N .v^ + X a X Ni v^) + 2A H X H v d v u + AmX Ni v 2 ^ = 0. 



free from the unwanted Nambu-Goldstone boson even if 
the singlet scalar S and/or iVj acquire VEV. For more 
details the reader is referred to Ref.[25, 43]. 

Minimization of the scalar potential (vide Eq.(ll)) 
leads to the following conditions: 



(14) 



(16) 

The solutions of the foregoing equation involve soft pa- 
rameters Am, Ah and m 2 . In fact these parameters can- 
not be much away from TeV values to have v s ~ TeV. 
In particular, the soft parameter Ah and Am are crucial 
to produce non zero VEV for the field S. Any consistent 
solution that yields v s ^ but . = requires \A\ > 3 
and also X H < 1, m > 100 GeV* m s > 100 GeV[25]. 
Similarly we also choose the couplings A s , AjVj in such a 
manner so that the condition for global minima is always 
satisfied. 



IV. NEUTRINO MASSES AND MIXING: 
R-PARITY CONSERVING NMSSM 

Let us now discuss in detail the generation of neutrino 
masses and mixing in this model. Note that this model is 
different from the models where MSSM is extended with 
three right-handed singlet neutrino superfields. This is 
because in those models the right handed neutrino mass 
scale is not tied up with the electroweak symmetry break- 
ing scale and is assumed to be very high (~ 10 15 GeV or 



I 

+ AM + v 2 ) + C ^^(v 2 v 2 u + ]>>£) 

i 

tan/3(-AifA 5 t>3 + + A h Xhv s )) + 2X H fl j v s vo i v^^ 

2v u (m Hu + X 2 H (v 2 + v 2 ) 9 -±±° 2 (v 2 -vl+Y, O 

i 

fl 3 fl k ^'N^N k + PvIu'v^vh + cotp{^X H X s v 2 + -XhXn.v 2 ^ + AhXhVs)) 
2A u fl j v j>i v$ j + 2fi i X Ni v s Vv j Vj i ., 

2v s (m 2 s + X 2 H {vj + v 2 ) + X s X H v d v u + X 2 N .v^ + -AmX s v a + ^ 2 s v 2 + -AgA^f^.) 
2A H XHVd.v u + AmXN z v 2 ^ + 2fjvv i v^. (AijUd + Ajv,-Uu), 

iv*M + - < + XX )) + 2A ^ v ^n d 

3 

2X H fJv d v sVfi . + 2fi k fi k vlvu, + 2fJX Nj v u VsV ii . + 2fJ ffv^v^v^, 
2vf l .(rn\f l , + AmX Nt v s + X Ni X H v d v u + -X Ni X s v 2 + X 2 N .v 2 + ^mAj^i^) 

yi i fi k vlvx k + 2A v fi i v u vz j +2X H fi i v d v s vz j +2X N J? 

I 



6 



so). 



Seesaw masses 



At the tree level, the (3 x 3) light neutrino mass matrix, 
that arises via the seesaw mechanism has a very well- 
known structure given by 



-m D M R 1 m 1 i D , 



where mo represents the lepton number conserving (3 x 
3) 'Dirac' mass matrix and Mr represents the lepton 
number violating (3 x 3) 'Majorana' mass matrix. Note 
that, after the EWSB, when the scalar component of S 
gets a VEV, in the effective Lagrangian we can assign a 
lepton number -1 for the fields N? and iV, (contained in 
the supcrfield Ni). The relevant part of the effective La- 
grangian which encompasses both neutrino and sneutrino 
fields is given by 



-L 



1 



eff = -{\mv s )NZNH + f» ]Vl v u NC + 



H.c. 



+ rnljiv* + + X 2 Ni v 2 )NN* 

+ {B^nNj + /OW; + B m NiNi + H.c), 



(18) 



where the coefficients have the following meaning 



fh i H — m z cos2(3, 



Bii = 
Bra = 



A v fijV u + X H f^v d v s , 

fijX Nj V u V s , 

^(X H X Ni v d v u + ^ s ^ Vs + Am\ Nl v s ). 



(19) 



It is easy to see from Eq.(18) that moij = fijV u and 
Mri = XniVs, which in turn provide neutrino masses at 
the tree level through Eq.(17). Note that, in Eq.(19) we 
have neglected a term ~ m 2 D in the expression for to?. 
since it is much smaller compared to the other terms. 

The tree level neutrino masses may receive dominant 
radiative corrections at the one-loop level. It has been 
shown in models of MSSM with right-handed neutrino 
supcrficlds, that the loop contributions can be as large as 
the tree level value, though the result depends on the soft 
SUSY breaking parameters [18, 20]. In P-parity conserv- 
ing scenarios the leading contribution to neutrino masses 
at the one-loop level arise from AL = 2 terms in the sneu- 
trino sector. These bilinear interaction terms involving 
the heavy right-handed sneutrinos fields Ni are B'^ViN*, 
and BRiNiNi as can be seen from Eq.(18). In associa- 
tion with the AL = term i.e., B^ViNj these AL = 2 
terms generate lepton number violating "Majorana" like 
mass terms (m 2 9 vC , + h.c.) for the left-handed sneutrinos. 
In fact, this can be seen as a scalar seesaw analogue of 



-X-.. 



Vi : 



> J J 



(17) FIG. 1: One-loop contribution to m„ when 



>0i = 0. 

Here Tn% i0 . represents the sneutrino "Majorana" mass term 
which generates the neutrino mass involving the sneutrino- 
neutralino loop. 



NjR(JI) 



'iR(iiy 



— >< < ■ -p. 

FIG. 2: The same one-loop contribution to m v as in Fig.l 
but represented in a different way. Here Njr and Nji are 
right handed sneutrino mass eigenstates which couple to Pi t j 
to produce the one-loop effective neutrino mass term. 



the usual fermionic seesaw mechanism to generate small 
masses for the light active neutrinos [20]. This effective 
Majorana sneutrino mass term in turn induces one-loop 
radiative corrections to neutrino Majorana masses via the 
self-energy diagram as shown in Fig.l. However, rather 
than computing the one-loop contribution to neutrino 
masses using the above method, we would choose a dif- 
ferent but more general procedure as explained below. 

We begin by decomposing the sneutrino fields in terms 
of real and imaginary components. Thus one has 



V2 



Ni = 



iR 



iNu 



V2 



(20) 



where the components UiR,NiR are the CP-even and 
Vu,Nii are the CP-odd scalar fields. The mass terms 
of these scalars may be evaluated using the definition 



M 



d 2 v 



d 2 v 



R,ij 



M 2 - = 



where $ represents a generic scalar field. Accordingly 
one obtains the following diagonal mass terms (assuming 
the right-chiral sneutrino states to be flavor diagonal) for 
the CP-even and CP-odd right-chiral sneutrinos: 



M 



R.NiNi 



m N,N* + X 



2 , 
Ni 1 
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+(^HXmVdV u + ^jv^^s + Am\ Nl v s ) 
-(XH^NiVdVu + ^NiKVs + Am\ Nl v s ). 



(22) 



Similarly, the interactions between N iR ^j^ and t'm^u) 
read as 



C 



li.v, V. — Jij 



fij^HVdVs + f"X Nj V u V s + A" />„, 



C P,u i N j = -fij^HVdVs + fij\NjV u V s - A v f"jV u . 

(23) 

The diagonal lcft-chiral sncutrino mass terms are shown 
in Eq.(19). As we can see, the off-diagonal terms involv- 
ing the lcft-chiral and right-chiral sncutrinos are much 



smaller compared to the diagonal terms since they are 
proportional to the small neutrino Yukawa couplings 
(/" ~ 10~ 6 ). Hence, we can compute the one-loop cor- 
rection to the neutrino mass due to the small mixing 
of the right-chiral sncutrinos with the left-chiral sneu- 
trinos. This is shown in Fig. 2. Note that the right- 
chiral sneutrino mass matrix contains bilinear terms like 
^H^NiVdVu, 2 iV " which are originated from the F- 
term contribution in the scalar potential. These are the 
new contributions to the right sncutrino masses in the 
present model and thus they are absent in seesaw mod- 
els of MSSM with only right handed neutrino supcrficlds. 
These terms will have important roles to play while cal- 
culating the one-loop correction to the neutrino mass 
matrix, even when the relevant soft breaking trilincar 
parameters are smaller. The loop contribution can be 
written as, 



J 



Mir 



fa 



X] TO x Q ( N a5 - tan0 w jV a4 ) 2 [ X CR.jCRjjhim^m^m^ , M^ JR ) 



,7=1,2,3 



,7=1,2,3 

where the integral I4 is given by, 

h{m„ iR , mo jR , m% a , m Xj ) 



d A q 



i(27r) 4 {q 2 - ml iR ){q 2 - m?^)^ 2 - m|J(g 2 - m 2 Xj ) 



(24) 



(25) 



Here Xj denotes right chiral sneutrino states Njr or Nj]. One can always evaluate I4 with the following analytical 
expressions, 



/ 4 (mi,m2,m3,m 4 ) 
I 3 (mi,m 2 ,m 3 ) 
J 2 (mi, m 2 ) 



1 



:[h(mi,m 2 ,m4) - J 3 (toi, m 2 , m 3 )], 
L 3 " "'4 

2^2(^1,7712) - /2(rni,m 3 )], 



1 



/Do 



(47r) 2 m{ - m\ 



lo S Tri- 



ms 



(26) 



Here, rn^ a represents the eigenvalues of the NMSSM neutralino mass matrix. In the weak interaction basis 
(jS, H^, H®, B°, W°^j , the mass matrix can be written as 



M = 



( Xsv s Xhv u X H v d \ 

Xhv u Xhv s -giVd/V2 g2Vd/V% 

^HVd X H v s giv u /V2 -giV n j\Fi 

-giVd/V2 giv u /V2 M ± 

V g 2 v d /V2 ~g2v u /V2 A/ 2 / 



(27) 



r 



The mixing matrix elements N a 5 and vV a 4 are the wino one-loop diagram (vide Fig. 2). Nevertheless, we would 
and bino component of the neutralino x a . The expres- consider a simplified scenario for illustration. In partic- 
sion (vide Eq.(24)) is the most general to compute the ular, we assume (i) identical values of Xi^i ( Ajy; = Xn) 
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for all three generations and (ii) soft-masses of the sneu- 
trinos (both v± and iVj) are flavor blind. This results into 
identical mass values for all three CP-even right chiral 
sneutrinos (Mft = Mjy ) and also for the three CP- 
odd states (Mftj = Mfj). With these assumptions, it 
is possible to factor out the flavor structure from Eq.(24) 
and denote the remaining as the loop factor (LF) which 
is merely a constant. Then the loop contribution can be 
cast into a convenient form given by 



(28) 



fe=i 



is related to Am? as m?~ = ^Am?[20]. The quantity /i 
is dchned as [i = \rv s . 

In order to reproduce the result in Eq.(34), we assumed 
that P/v ) to jvjv* ^ Mr and A v > Pjv. Now, in addition 
if we assume Mr > A v , the last term becomes negligible 
compared to the other terms in the expression Eq.(35) 
and this keeps only the terms to leading order in 1/Mr. 
However, this is not always true as all soft SUSY breaking 
mass parameters as well as the right handed neutrino 
masses may have similar magnitudes as in the present 
scenario. Hence, rather than using Eq.(34), we evaluate 
the neutrino mass terms corrected up to one loop order, 
from 



where 



LF = xX] m *<*(^5-tan0 w iV a 4) 2 



(I 4 (mp,mp,m^ a , M^ r )C% - h(mp,mp,m% a , M^^Cp). 

(29) 

Here Cr and Cp represent the coefficients of /y in Eq. 
(23) and given as 

Cr = X H v d v s + X N v u v s + A v v u , (30) 
Cp = -XHVdV s + \nv u v s — A v v u . (31) 

Let us note that the coefficient Bjh can be written as 

B Rl = B N M R (32) 

where 



Pjv 

Mr = \ N V. 



^(X H VdV u /v s + ^y^- + Am), 



(33) 



Consequently the one-loop contribution can be cast 
into the well known form[18, 20] 



{loop) 



f(Va) 



g%Ampij 
327r 2 cos 2 8w 

(1-jm) 2 



(34) 



where y a = m?/m| and N a ^ = N a5 cos 6w — N ai sin. 6yy 
is the neutralino mixing matrix element and to order in 
l/Mft the left sneutrino mass difference relative to the 
light neutrino mass is given by 



2(A V + n cot f3 - B N - 



B N (A„+fi cot/3)' 
Ml 



) 



(35) 

Here we have used the relation Am? = 2mpAmp and rap 
is an average left-sneutrino mass. In the present case all 
left handed sneutrino soft masses are assumed to be iden- 
tical. The sneutrino Majorana mass m?p shown in Fig.l 



Mr 



k=l 



(36) 



Clearly, the coefficient of the loop contribution shifts the 
tree level neutrino masses by a constant amount. This 
coefficient involves the soft SUSY breaking parameters 
and in this work we explore the effect of these parameters 
on the neutrino mass matrix. 

This simple structure of the neutrino mass matrix 
(vide Eq.(36)) can indeed be very helpful to examine 
the neutrino mixing pattern. In particular, we are in- 
terested to explore the conditions which could yield the 
mixing matrix into a tri-bimaximal structure. Thus we 
compare Eq.(36) with Eq.(4), where the latter provides 
with the neutrino mass matrix consistent with the tri- 
bimaximal mixing pattern. Then, with a symmetric neu- 
trino Yukawa matrix, neutrino masses can be evaluated 
using the following expressions: 



-mi + -m 2 
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(m 1+ 2m 2 + 3m 3 ) = C[(ft 2 ) 2 + (ft 2 ) 2 + (f% 3 )% 

i(-m 1 + m 2 ) = C[/ri/r 2 + AV22 + AV23L 
h mi -m 2 ) = CIAVra + AVas + AVas], 



-(-m a - 2m 2 + 3m 3 ) = C[AVi3 + A2A3 + /M- 

(37) 

_ 2 

Here the constant C is defined as (-^+LF). As a simple 
choice we consider, f 22 = fm an d a ^ so A2 — A3 = to 
obtain the solutions. This choice, coupled with the con- 
sistency condition Ai = AT2 — A?3> leads to the following 
solutions of the neutrino spectra 



mi 



m 3 



m 2 = (—± + LF)(K 1 f, 
Mr 



v \2 



Mi 



IP)(2/ 2 Wn) 



(38) 
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It is obvious that the mass pattern as depicted above sat- 
isfies the desired tri-bimaximal structure of the neutrino 
mixing. The mass terms as expected, contain tree level 
contributions which are always negative. On the other 
hand, the loop contribution can go both ways depending 
on the sign of the soft SUSY breaking parameters. For a 
large Br, which primarily depends on Am, the radiative 
correction to the neutrino masses could be enhanced to 
supersede the tree level results[18]. 

Before presenting the numerical results a few com- 
ments regarding the lepton flavor violating (LFV) pro- 
cesses are in order. Recall that we assume flavor diag- 
onal mass terms for the left and right chiral sneutrinos. 
The loop induced processes like fx —> &y, r — > or 
t — > fX'y can get contributions primarily via the couplings 
B"j or B% (see Eq.(18) and (19)). Clearly, any such con- 
tribution at the leading order would involve a product 
of two small neutrino Yukawa couplings ffj and are ex- 
pected to be very suppressed. Moreover, our assumption 
/12 = fia = would lead to vanishing contributions for 
the processes fj, — > ej and r — > e^j in this model. 

We now explore whether the obtained mass pattern 
could fit with the different hierarchical structure that we 
know so far. In particular, we show our numerical results 
to identify the regions in the parameter space consistent 
with the normal, inverted and quasi-degenerate neutrino 
mass pattern. In the numerical computation we choose 
different soft parameters and couplings in such a way, 
that the proper minima condition of the scalar potential 
is always satisfied [25]. 

The choices of various parameters are listed below. 
The value of tan/3 is taken to be equal to 10. In addition 
to that, other parameter choices are 

(I) Superpotential parameters: A^ = —0.3, A s = 0.6, 
Ajvi = ^N2 = Ajv3 = A;v = 0.2, 

and 

(II) Soft SUSY breaking parameters: m s = 100 GeV, 
mfj.ff, = 300 GeV, m 9i = 100 GeV, A H = 100 GeV, 
A v = 1000 GeV. 

Apart from the above parameters which are fixed to 
the quoted values, we have also varied the parameter 
Am in the calculation. This would cause changes in v s 
(vide Eq.(16)), which in turn produces variation in the 
neutrino spectrum. We list the values of Am and v s in 
table I. 



Am (GeV) 


-600.0 


-800.0 


-1000.0 


-1200.0 


v s (GeV) 


927.56 


1280.76 


1625.18 


1965.67 



TABLE I: Different values of v s corresponding to the different 
values of the coupling parameter Am. 



B. Different Neutrino Spectra: 

The two mass-squared differences shown in Eq.(l) indi- 
cate three possible neutrino mass hierarchies [44] , namely 

1. Normal Hierarchy: this neutrino mass pattern can 
be established if mi,m 2 and 7 77,3 are r elated with 
the observables y'Am^ and VdAm|^[ as 



mi ps m 2 ~ y Amli, m 3 ~ 1 Am§ 2 | . (39) 

However, in principle m\ can also be much smaller 
than m 2 or even be zero. Since in this case 7713 is 
much greater than both mi and mi, we can ap- 
proximately use the relation shown in Eq.(39) for 
illustration. 

2. Inverted Hierarchy: this hierarchical scenario can 
be achieved if one chooses 



mi ps m 2 ~ Y |Amg 2 |, m 3 <Y|Am§ 2 |. (40) 

We assume the maximum possible value for 7713 to 
be ~ 0.01 cV while the minimum value could be 
vanishing. Obviously, the solar mass squared dif- 
ference Am^i will come from the small mass split- 
ting between m 2 and mi, where Am^j <C m 2 ,mi. 
Hence, for a simple minded analysis we can assume 
that m 2 = m\. 

3. Degenerate Masses: finally this scenario is defined 

by 



mi ps m 2 ps m 3 > J lAm^l, (41) 

Here we assume that the upper bound of the neu- 
trino masses could be 0.33 cV, which comes from 
the cosmological observations. The lower bound is 
chosen to be 0.1 eV. 

In Fig3, three neutrino mass eigenvalues mi,m 2 ,m 3 , 
consistent with the normal hierarchical pattern, arc plot- 
ted as functions of neutrino Yukawa couplings. The dif- 
ference in the contours manifests how the neutrino masses 
depend on the soft bilinear coupling parameter (Am). 
The variation occurs, as v s depends on (Am), thereby 
acquiring a different value at the global minima which 
has already been mentioned in Table I. In particular \v s \ 
always increases as we increase \Am\ parameter which in 
turn increases the right handed neutrino masses. This 
results into a smaller value for m l J ee . On the other 
hand, loop correction does not increase appreciably by 
this small variation of Am if A v is around TeV scale as 
we will discuss later. We should note here that neutrino 
loop correction is always an order of magnitude smaller 
compared to the tree level value for the parameters we 
have chosen. Thus with increase in Am parameter, one 
requires large values of Yukawa couplings to satisfy the 
neutrino data. The red zone in each contour (vide Fig3) 
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FIG. 3: Normal hierarchy: variation of m„ with the Yukawa 
parameters. The red (solid) segments denote the range of the 
Yukawa parameters that satisfy neutrino data. Each contour 
represents a separate set of v a and Am, as given in Table. I. 
All mass parameters are in GeV. 



FIG. 4: Variation of 777,3 with the corresponding Yukawa pa- 
rameter is shown for (a) inverted hierarchical mass pattern 
and also for (b) the degenerate spectrum. All mass parame- 
ters are in GeV. 



represents the range of the Yukawa couplings that can 
satisfy the neutrino data. 

In case of inverted hierarchy, we have shown the vari- 
ation of 77z 3 with the respective Yukawa couplings in 
Fig. 4(a). The other mass parameters mi,m2 depend on 
the Yukawa coupling f" x , but that can be estimated from 
the Fig. 3(b) if in that plot we replace wis in the y-axis by 
m\/m2 and 2f 22 — fii in the x-axis by [vide Eq.(38)). 
In fact knowing the value of the Yukawa coupling 
would allow us to determine the coupling f 22 . 



The Fig. 4(b) depicts the variation of 777,3 with the 
Yukawa coupling f 22 for quasi-degenerate mass scenario. 
In this scenario, the neutrino spectrum is approximately 
degenerate i.e., mi, 7712 and 777,3 turn out to be almost 
identical if one chooses f 2i much smaller compared to the 
diagonal Yukawa coupling f 22 , which essentially means 
that 

Finally a few comments on the dependence of the one- 
loop contribution to the neutrino mass on the soft SUSY 
breaking parameters Am and A v . The loop contribution 
is always suppressed unless the parameter A v is sum- 
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1000 3000 5000 7000 9000 11000 13000 15000 

A v (GeV) 

FIG. 5: Variation of the m l °° p with A" in the normal hierar- 
chical scenario. All mass parameters are in GeV. 



ciently large as can be seen from Fig. 5. As for illustration, 
the Yukawa couplings are chosen as /fj = 1.75 x 10 -7 
and f% 2 = /33 = 3-95 x 10~ 7 . Similarly, we choose Mi 
= 60 GeV and M 2 = 120 GeV, where M x and M 2 are 
the U(l) and SU(2) gaugino mass parameters, respec- 
tively. For larger values of electroweak gaugino masses, 
the one-loop contribution would be reduced further. For 
Am = -1 TcV, higher A v values ~ 13 TcV can satisfy the 
current neutrino data. However, even if A v is ~ 13 TcV, 
the quantity Ayj u is very small, i.e., ~ 10~ 2 GeV. Note 
that for such a choice of the parameter space, the tree 
level values of the neutrino masses are not sufficient to 
accommodate the three flavor global neutrino data. In- 
creasing the value of Am requires relatively smaller value 
of A v (~ 7 TeV) to reproduce the neutrino data. It is 
very important to point out that, for a fixed A s and X^i 
one cannot increase the trilinear coupling parameter Am 
to an arbitrary high value as the right-chiral sncutrinos 
may turn out to be tachyonic. Thus, a relatively larger 
soft trilinear parameter A v is required to enhance the 
one-loop contribution to neutrino masses. The require- 
ment of a large A v can be understood from the following 
discussion. 

• The one-loop contribution to the neutrino mass 
originating from the mass splitting in the left- 
handed sneutrinos depends on the parameters /x, 
A v and Bn as can be seen from Eqs.(34) and 
(35). It has been argued in Ref.[18], that in order 
to have the one-loop contribution to the neutrino 
mass comparable to its tree level value, the ratio 
Am^./m,, should be ~ 10 3 . 

• Substituting fi = Xhv s and the expression for Bn 
from Eq.(33), we may write Am vi j /m„ ~ 2(A V + 



X H v s cot p - (\X s v s + Am/2 + X H v d v u /2v s )(l + 
(Ay + X H v s cot P) 2 )/Mft)/my. We can see from 
the above expression that one may increase either 
Am or A v parameter to enhance the one-loop con- 
tribution to make it countable. But in the present 
context, raising the soft parameter Am alone would 
not serve the purpose. This is because the VEV v s 
increases significantly with \Am\ (vide Table. I). 
Thus there is always a partial cancellation between 
different terms in the above expression for the left 
sneutrino mass splitting. In particular, the effective 
bilinear coupling Bn is reduced because of this par- 
tial cancellation. In addition, we choose the sign of 
the coupling Xh as negative in order to determine 
the correct global minima. This also causes a par- 
tial cancellation between various terms, but to a 
lesser extent. Considering this cancellation effect in 
mind, it is easy to check that the ratio Ammj/my 
always reside near the value ~ 10 with the soft pa- 
rameters A v and Bn around the TeV scale. 

• Now, as mentioned above, the trilinear coupling 
parameter Am is restricted if one does not want 
the right chiral sncutrinos to become tachyonic. 
Of course this depends on the choice of the soft 
"Dirac" mass term mfjfj, of the Ns, which we have 
chosen to have a quite moderate value (300 GeV) 
in this case. However, the parameter A v can be 
pushed to a reasonably high value without affect- 
ing any other results. This explains why a large A v 
parameter is required to make the one-loop contri- 
bution to the neutrino mass comparable to its tree 
level value. 



V. SIGNATURES AT LHC 

It is extremely important to investigate the possible 
signatures of this TcV scale seesaw mechanism at the 
LHC. One of the search strategics could be to produce 
the right-handed neutrino N (or the corresponding right- 
handed sneutrino N) with a large enough cross-section 
and then look at the decay branching ratios in differ- 
ent available modes. However, in this type of models the 
production of TeV scale right-handed neutrinos (or sneu- 
trinos) at the LHC is suppressed 1 by the light neutrino 
mass [46]. Nevertheless, it is possible to construct mod- 
els where the production mechanism of the right-handed 
neutrino (sneutrino) can be decoupled from the neutrino 
mass generation. For example, extended gauge symme- 
tries such as U(1)b-l or SU(2)r may offer extra gauge 
bosons near the TeV scale whose couplings to quarks and 
the right-handed neutrino (sneutrino) are unsuppressed 



A very recent analysis along with the discovery potential at the 
LHC is presented in Ref. [45]. 
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[47]. In such models a single or a pair of right-handed 
neutrinos can be produced with large cross sections lead- 
ing to dilepton signals (same-sign) with no missing en- 
ergy (see the first reference of [17] and [48-51]), trilcpton 
signals [52] or four-lcpton signals [53-55]. 

In the context of the present model the left-sneutrino 
"Major ana" mass term can lead to oscillation between 
the left-chiral sneutrino and the corresponding anti- 
sneutrino [18, 56, 57]. This can be interpreted as the ob- 
servation of a sneutrino decaying into a final-state with a 
"wrong-sign" charged lepton. In order to have a large os- 
cillation probability the total decay width T of the sneu- 
trino/antisncutrino and the mass splitting Am must be 
of the same order. Since Am is constrained by the neu- 
trino data, one needs a very small total decay width of the 
sneutrino/antisneutrino. It has been shown in [18] that 
this can be achieved in a scenario where the lighter stau 
is long-lived and the left-chiral sneutrino can only have 
3-body decay modes involving the lighter stau in the final 
states. This can lead to signals such as like-sign dilep- 
tons, single charged lepton plus like-sign di-staus (lead- 
ing to heavily ionizing charged tracks) or like-sign di-stau 
charged tracks at future linear colliders [18, 58, 59] or at 
the LHC [60]. The resulting charge asymmetry of the 
final states can be measured to get an estimate of the 
sneutrino- antisncutrino oscillation probability [60]. In 
addition, for a very small sneutrino decay width one can 
also observe a displaced vertex in the detector. However, 
a detailed study of such signals in the context of the 
present model is beyond the scope of the present paper. 

In comparison, now we discuss briefly the signatures 
of R-parity violating models in general. In models with 
spontaneous violation of R-parity, the singlet sneutrino 
vacuum expectation value leads to the existence of a Ma- 
joron which is an additional source of missing energy. 
This can change the decay pattern of the lightest Higgs 
and the lightest neutralino with the corresponding sig- 
natures at the LHC. For more details and the relevant 
references the reader is referred to Ref.[46]. In the case 
of bilinear R-parity violation, the ratios of certain de- 
cay branching ratios of the LSP show very nice correla- 
tion with the neutrino mixing angles. This can lead to 
very interesting signatures at the LHC where compara- 
ble numbers of events with muons and taus, respectively, 
can be observed in the final state [9-15]. 

From the above discussion we see that the canonical 
type-I supersymmetric seesaw case that we have consid- 
ered in this paper has characteristic signatures which can 
be tested at the LHC. At the same time one can also dis- 
tinguish the predictions of this model with those of the 
models with spontaneous or bilinear R-parity violating 
scenarios. 



VI. CONCLUSIONS 

We have studied the neutrino masses and mixing in 
an R-parity conserving supersymmetric standard model 



with three right handed neutrino superfields iV, and an- 
other gauge singlet superfield S. This model is similar 
to the next-to-minimal supersymmetric standard model 
(NMSSM), where the scalar component of S gets a VEV 
to generate a /i-term of correct order of magnitude. In ad- 
dition, the same VEV also generates TeV scale Majorana 
masses for the right handed neutrinos. The small neu- 
trino masses are generated at the tree level by the usual 
seesaw mechanism at the TeV scale. We also calculate 
the one-loop contribution to the neutrino mass matrix 
and investigate the constraints on the model parameters 
to produce the tri-bimaximal pattern of neutrino mixing 
for three different neutrino mass hierarchies. Neutrino 
mass matrix gets contribution at the one-loop level con- 
trolled by the sneutrino "Majorana" mass terms. We 
show that the one-loop contribution can be important 
for certain choices of the soft SUSY breaking parame- 
ters. This we have demonstrated by evaluating the one- 
loop contribution in two different ways. In particular, 
we observe that the one-loop contributions can be signif- 
icant when the soft SUSY breaking trilinear parameter 
A v f v is - O(10- 3 GeV) with A v ~ 10 TeV. This ob- 
servation is quite robust and does not change much if 
one introduces a small #13 in the neutrino sector. Our 
choice of neutrino Yukawa couplings also predict vanish- 
ing contributions to the lepton flavor violating processes 
fi — > ej and r — > ej as well as an extremely suppressed 
contribution to r — > fij. 

As has been stated earlier, it is also possible to have 
non-zero vacuum expectation values for the left and right 
chiral sneutrinos. In that case, R-parity is violated spon- 
taneously. The neutrino mass matrix can have contri- 
butions from two different sources, namely, the effective 
bilinear R-parity violating interactions and the TeV-scale 
seesaw mechanism. One-loop contributions to the neu- 
trino mass matrix can be very important in this case 
too. However, the tree level and one-loop calculations 
are rather involved and require a separate discussion al- 
together. We plan to present these results in a subsequent 
paper [42]. 

The characteristic signatures of this model at the 
LHC include like-sign dilepton (without missing energy), 
trilcpton or four lepton final states as well as single lep- 
ton plus two heavily ionizing charged tracks or only two 
heavily ionizing charged tracks stemming from long-lived 
staus. By looking at these signals one can possibly dis- 
tinguish this model from the models of spontaneous or 
bilinear R-parity violation. 
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